Determination of the Axial- Vector Weak Coupling Constant with Ultracold Neutrons 
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A precise measurement of the neutron decay /3-asymmetry Ao has been carried out using polar- 
ized ultracold neutrons (UCN) from the pulsed spallation UCN source at the Los Alamos Neu- 
tron Science Center (LANSCE). Combining data obtained in 2008 and 2009, we report Aq — 
—0.11966 ± 0.00089t 000 1 40 , from which we determine the ratio of the axial-vector to vector weak 



coupling of the nucleon qa/ 



-1.27590 



+0.00409 
0.00445- 



PACS numbers: 14.20.Dh, 12.15.Ff, 12.15.Hh, 23.40.Bw 



The axial- vector weak coupling constant, gA, plays an 
important role in our understanding of the nucleon spin 
and flavor structure |l|, [2J . It is a central target for high 



Fermi coupling constant Gf, measured to 5 ppm using 



and an essential 
gA is also im- 



precision lattice QCD calculations [3, 

parameter in effective field theories 

portant in a variety of astrophysical processes, including 

solar fusion reaction rates [6J. 

The angular distribution of emitted electrons from po- 
larized neutron decay can be expressed as W(E) oc 1 + 
- (P) A(E) cos 8 , where A{E) specifies the /3-asymmctry 
versus electron energy E, v is the electron velocity, c is 
the speed of light, (P) is the mean polarization, and 8 
is the angle between the neutron spin and the electron 
emission direction [7[. The leading order value of A(E), 
A , is given by 



muon decay [Uj, the CKM matrix element |V^d|, mea- 
sured to 225 ppm using + — > + supcrallowed decays 
[8j, and the neutron lifetime, measured to 0.9% [12J [131 ] . 
Thus, a measurement of the /3-asymmetry permits direct 
determination of gA , as well as a robust test of the con- 
sistency of measured neutron /3-decay observablcs under 
the Standard Model. 

In order to obtain A(E), one must determine the po- 
larization of the neutron beam and control all sources 
of systematic uncertainty due to backgrounds, including 
those produced by the neutrons themselves, the detec- 
tor response, and electron-event reconstruction. All pre- 
vious precise measurements of the /3-asymmetry 
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-2(A 2 - |A|) 
1 + 3A 2 ' 



(1) 



where A = gA/ ' gv and gv is the vector weak coupling 
constant with gv = 1 under the conserved vector current 
(CVC) hypothesis of the Standard Model 0. Higher 
order terms in A{E) are at the 1% level, and can be 
calculated precisely under the Standard Model [9j, llOj . 
gA can also be indirectly determined by combining the 



have been performed with cold-neutron beams and have 
shown a range of results much wider than the reported 
uncertainties [12j. Our measurement, UCNA, utilizes ul- 
tracold neutrons (neutrons with kinetic energy less than 
200 neV) and controls key systematic uncertainties: neu- 
tron polarization and neutron-generated backgrounds. In 
2007, we carried out a proof-of-principle /3-asymmetry 
measurement |lSJ] . At present, the UCNA experiment 
is characterized by neutron polarizations greater than 
99.48% and neutron-generated backgrounds that produce 



corrections to the asymmetry below the 0.02% level. 

Some of the experimental details of UCNA are ex- 
plained in [3. We used the UCN source at the LANSCE 
accelerator at Los Alamos National Laboratory [19( . The 
UCN were polarized by a 7 T primary polarizer coupled 
to an adiabatic fast passage (AFP) spin flipper to control 
the spin state lg, [20j • Polarized UCN entered the super- 
conducting spectrometer (SCS) [21| and were confined in 
a 3 m long, 12.4 cm diameter electropolished Cu tube 
(decay trap) with variable thickness mylar endcaps. The 
inside surface of each endcap was coated with 200 nm 
of Be. A 1 T magnetic field was oriented parallel to the 
decay trap, along which decay electrons spiraled toward 
one of two identical electron detector packages, each cov- 
ering a 2-7T sr angular hemisphere. Each detector package 
consisted of a low-pressure multiwire proportional cham- 
ber (MWPC) [22| backed by a plastic scintillator, with 
scintillation light measured by four photomultiplicr tubes 
(PMT). Each MWPC had thin front and back mylar win- 
dows which separated low-Z chamber gas (neopentane) 
from the spectrometer vacuum. To study key systemat- 
ics due to electron energy loss and backscattering in the 
windows, we operated the experiment in four different ge- 
ometries with different decay trap endcaps and MWPC 
window thicknesses, A: 0.7 and 25, B: 13.2 and 25, C: 0.7 
and 6, and D [29|: 0.7 and 6 fim, respectively. 

Cosmic-ray muon backgrounds were identified by a 
combination of plastic scintillator veto paddles and sealed 
drift tube assemblies [23j surrounding the electron detec- 
tors. 7-ray backgrounds were vetoed by a coincidence 
between the MWPC and the main /3-scintillator. 

A gate valve separated the UCN source from the ex- 
perimental apparatus. A typical run unit consisted of a 
background run (gate valve closed), a /3-dccay run (gate 
valve open), and a UCN depolarization run to measure 
the equilibrium UCN polarization for the accompanying 
/3-decay run. The UCN spins were flipped back and forth 
(while the magnetic field in the SCS was held fixed) 
between run units (Fig. [TJ Panel (a)), which partially 
canceled systematic rate variation over the period of a 
spin cycle (~1.5 hour). During a depolarization run, the 
guide serving as input to the 7 T polarizing field was first 
connected to a UCN detector [2J] so that UCN exiting 
the experiment could be counted, while the gate valve 
was closed and the proton pulses [19| were discontinued. 
This cleaning phase, which lasted 25 s, produced a sig- 
nal in the UCN detector proportional to the number of 
correctly polarized UCN present in the experimental ge- 
ometry at the end of the /3-decay measurement interval. 
Following the cleaning phase, the state of the spin flipper 
was changed, allowing only incorrectly polarized UCN re- 
maining downstream of the spin flipper to pass through 
the 7 T polarizing field and be counted. Counting during 
this unloading phase was performed for ^200 s in order 
to measure background as well as incorrectly polarized 
UCNs. Since the measured depolarization was consistent 



with zero at the la level, we folded together statistical 
and systematic errors to produce a global polarization 
lower limit of 99.48% at the 68% CL [U, covering all 
four geometries and both polarization states. 

The experimental triggers were formed by requiring 
at least 2-of-4 PMT signals over threshold in either of 
the scintillator detectors. Electron positions were deter- 
mined with the MWPC to an accuracy of better than 2 
mm based on the distribution of charge on two perpen- 
dicular cathode grids in the MWPC [221 ■ A fiducial cut 
of r < 45 mm was placed on the trigger side to reduce 
background and to eliminate electrons that could strike 
the decay trap walls. 

Reconstructed event energies E rocon were measured us- 
ing the signals from the scintillator PMTs as calibrated 
with conversion electron sources ( 109 Cd, 139 Ce, 113 Sn, 
85 Sr, and 207 Bi). The position-dependence of the re- 
sponse of each PMT was mapped out by comparing the 
neutron beta decay spectrum endpoint observed at dif- 
ferent positions. The energy reconstruction uncertainty 
was determined to be the larger of ±5 keV or ±2.5%, 
which covered the uncertainty in the position response of 
the scintillator, as well as possible variation of the energy 
response allowed by the calibration data. 

The PMT gains were monitored based on frequent cali- 
brations with a removable 113 Sn calibration source, which 
also measured the energy resolution of the system (~ 400 
photoelectrons per MeV), and by observing shifts in the 
minimum-ionizing peak of cosmic-ray muons during /3- 
decay and background runs. 

The majority of the /3-decay events were single detec- 
tor triggers. However, due to electron backscattering, a 
small fraction of the events, varied between 1.7% and 
3.4% for the four geometries, triggered both scintilla- 
tors and another small fraction (~2%) were detected by 
both MWPCs, but triggered only one of the scintilla- 
tors. In the first case, the initial direction of the electron 
could be determined by the relative timing of the trig- 
gers, while in the second case a fixed cut (4.1 keV) or a 
likelihood function based on the energy loss in the trigger 
side MWPC yielded an identification efficiency of ~80% 
based on Monte Carlo calculations (discussed later). 

In addition to the ambient backgrounds (measured 
with the UCN gate valve closed and suppressed by the 
pulsed nature of the UCN source [19( ) , which were sub- 
tracted run-by-run, neutron captures in the vicinity of 
the detectors could create prompt 7's with energies up 
to ~ 8 MeV, generating an irreducible background in 
the experiment. This background was significantly sup- 
pressed, compared to cold neutron beam experiments, by 
the relatively low density and low capture and upscatter 
probability of neutrons in and around the spectrometer. 
Combining direct measurements with Monte Carlo cal- 
culations, we obtained an upper limit of 0.02% on the 
correction to the asymmetry. 

For each run, events were sorted into 25 keV E lecon 



bins from to 1200 keV and assigned an initial direc- 
tion. The rates in the two detectors were then computed 
based on the experiment live time. We applied sepa- 
rate spin-dependent blinding factors to the two detector 
rates, effectively adding an unknown scaling factor to the 
measured asymmetry that was constrained to be within 
1.00 ± 0.05. After determination of all cuts, corrections 
and uncertainties, this factor was removed. For each 
/3-decay/background run pair, the background rate was 
subtracted from the /3-decay-run rate bin by bin. The 
reconstructed energy spectrum (background subtracted, 
averaged over the two spin states) is shown in Panel (b) of 
Fig.Q] overlaid with the measured background. The S/B 
is about 40 in our analysis energy window between 275 
and 625 keV (discussed later). Also overlaid is the Monte 
Carlo-predicted reconstructed energy spectrum, with all 
detector effects (efficiencies, resolutions, etc.) taken into 
account. The systematic effect due to the small discrep- 
ancy between the two spectra is well covered by the en- 
ergy reconstruction uncertainty in Table |U 

For a given geometry, a "super-ratio" of count rates 
among the two detectors and UCN spin states was cal- 
culated (as defined in [la]), from which the raw mea- 
sured asymmetry was determined (see also Fig. [1] Panel 
(a)). To extract Aq, we first multiplied the raw measured 
asymmetry by l/(v/c) in each energ y b in to remove the 
strongest energy dependence. As in [la], two scattering- 
related effects dominated subsequent systematic correc- 
tions: the residual backscattcring correction and the an- 
gle effect. In addition to a small residual correction due 
to incorrect identification of the initial electron direc- 
tion for the measured electron backscatters (where both 
detectors observed the electron), there were corrections 
for backscattcring from the decay trap windows and the 
front windows of the MWPC that could not be identified 
experimentally Angle effects arose from the fact that 
the energy loss of an electron in the thin windows was 
strongly angle-dependent. Low-energy, large pitch angle 
electrons were more likely to fall below the scintillator 
threshold, leading to a suppression of the acceptance at 
large angles ((cos 8) deviating from 1/2). Both of these 
effects were evaluated with two independent simulation 
programs: Penelope [H] and GEANT4 [2(|. The re- 
sulting corrections for all four geometries are shown in 
Table HI Based on the observed difference between the 
calculations and the data, we assigned an uncertainty of 
30% to the backscattering correction and 25% to the an- 
gle effect correction. 

Recoil-order corrections to A(E) (see also la]) were 
calculated within the context of the Standard Model ac- 
cording to the formalism of [9j, [lOJ , leading to a correction 
of —1.79 ± 0.03% to Aq. The value for the radiative cor- 
rection to Aq was taken from the calculations of [271 ] . 
yielding a small theoretical correction of 0.10 ± 0.05%. 

Applying all corrections mentioned above, the ex- 
tracted Aq is plotted against E rccon (all geometries com- 
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FIG. 1: (Color online) Panel (a): schematic of the experiment 
and definition of the asymmetry. Panel (b) : background sub- 
tracted electron -E re con spectrum (solid circles with the uncer- 
tainty of -Erccon reflected by the horizontal error bars), com- 
bining both sides of the detectors and two spin states, overlaid 
with the Monte Carlo spectrum (histogram). The open cir- 
cles represent the measured ambient background spectrum. 
Panel (c): Ao vs. -E r0 con, combining all four geometries. The 
horizontal line represents the extracted Aq within [275, 625] 
keV. Drawn on the same scale below the graph are four sets 
of bands representing the sum of energy-dependent backscat- 
tering and angle effect corrections (light color), positive sign 
indicating a larger \Aq\, and their uncertainties (dark) for the 
four geometries. The positive (negative) correction at low 
(high) energy is a consequence of the backscattering (angle 
effect) dominating in this energy region. The number of j3 
events for each geometry is also indicated in the figure. 



bined) in Panel (c) of Fig [T] Energy-dependent correc- 
tions (backscattering and angle effects) and their uncer- 
tainty are indicated as bands in the figures. The final Aq 
is obtained from a constant fit over a range of energy [28( . 
The energy window, 275 to 625 keV, was chosen to op- 
timize combined statistical and systematic uncertainties 
before unblinding the asymmetries. The value of Aq was 
insensitive to the choice of energy window, with the vari- 
ation less than 15% of the statistical uncertainty for win- 
dows between 150 and 750 keV. 

The experimental uncertainties and systematic cor- 
rections to Aq arc summarized in Table [TJ Geometry- 
dependent systematic uncertainties (backscattering, an- 



TABLE I: Summary of experimental corrections and un- 
certainties in % (all fractional to Aq). Upper: geometry- 
independent effects. Lower: geometry-dependent effects (first 
value=correction, second value=uncertainty in each column), 
with (Tstat, Aback, A ang and emwpc referring to statistical un- 
certainty, backscattering correction, angle effect correction, 
and the uncertainty associated with MWPC inefficiency, re- 
spectively. 



Geometry-independent effect corr. (%) 

Polarization 

Field non-uniformity 

Rate dependent gain shift 

Gain fluctuation 

Deadtime 

Energy reconstruction 

UCN-induced background 

Muon veto efficiency 

Live time uncertainty 

Fiducial cut 



unc. (%) 

TU752 

-0 

+0.20 

-0 

0.08 
0.20 
0.01 
0.47 
0.02 
0.30 
0.24 
0.24 





Geometry-dependent effect 






A(%) 


B(%) 


C (%) 


D(%) 


Cstat 


n/a 1.23 


n/a 1.22 


n/a 2.00 


n/a 2.10 


Aback 


1.34 ±0.40 


4.32 ±1.30 


1.07 ±0.32 


1.08 ±0.32 


^ang 


-1.81 ±0.45 


-3.22 ±0.81 


-0.60 ±0.15 


-0.36 ±0.09 


EMWPC 


0.02 


0.01 


0.16 


0.5 



gle effect, and MWPC inefficiency) arc treated as com- 
pletely correlated among the different geometries. Com- 
bining the four geometries, we find Aq = —0.11966 ± 



0.00089 



+0.00123 
-0.00140 



(with a x l v °f 2.4/3), where the 
first uncertainty is statistical and the second system- 
atic [28|. Based on Eq. ([T]), we then determine gA/gv = 
- gA, where the second equality as- 



1.27590_ 00445 
sumes CVC|8j. 




This worl^ (2010) 
[14] (2002) 
[15] (1997) 
[16] (1997) 
[17] (1986) 



-0.124 -0.122 -0.120 -0.118 -0.116 -0.114 -0.112 -0.110 -0.108 -0.106 -0.104 
A. 



FIG. 2: (Color online) Ideogram comparing the previous mea- 
surements that are included in the 2009 Particle Data Group 
best value for Aq along with the new result reported here. 
The dashed line is the ideogram with the previous four data 
points, and the solid line shows the result of including the 
present data. 

Our result for Aq is compared with the world data |l4J - 
17] in Fig. [5J Our result is in good agreement with the 



pendent of the CKM matrix element \V u d\- This strongly 
motivates UCN-based measurements of the neutron /3- 
asymmetry, where the key neutron-related systematic un- 
certainties can be reduced below the 0.1% level. Recently 
demonstrated improvements to the UCN source and re- 
finement of the energy response and gain monitoring will 
permit the collection of a much larger data set and the 
reduction of all major systematic uncertainties. 
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Energy Office of Nuclear Physics, National Science Foun- 
dation (NSF-0555674, NSF-0855538, NSF-0653222), and 
the Los Alamos National Laboratory LDRD program. 
We gratefully acknowledge the support of LANSCE and 
AOT divisions of Los Alamos National Lab. 



most recent and precise result for Aq [L4j. We note that 
the direct extraction of gA/gv from the /^-asymmetry is, 
unlike extraction from the neutron lifetime [12:1 301, inde- 
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